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SUMMARY 


f^R^ults  of  a  numerical  study  of  the  natural  frequencies  and  modes  of 
vibration  of  the  pressure  prestressed  toroidal  shell  are  presented.  The 
analysis  is  based  upon  a  linearized  theory  of  vibrations  of  prestressed 
shells.  The  frequencies  and  mode  shapes  are  obtained  j^y^rial  and 
error  in  the  Ilolzer  fashion.  The  effects  of  wall  bending  stiffness  on  the 
frequencies  of  shells  under  varying  degrees  of  prestress  are  shown,  jj  ^ 
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Geometry  and  Notation 

Axisymmetric  Modes,  Symmetric,  k  =  0.002 
Axisymmetric  Modes,  Antisymmetric,  k  =  0.002 
n  =  1  Modes,  Symmetric,  k  =  0.002 
n  =  1  Modes,  Antisymmetric,  *  =  0. 002 
n  =  2  Modes,  Symmetric,  «■  =  0.002 
n  =  2  Modes,  Antisymmetric,  k  =  0.002 

Effect  of  Bending  Stiffness  on  Mode  Shapes,  n  =  0,  «•  =  0.  002,  e  =  0.  75 
Axisymmetric  Modes,  Symmetric,  «•  =  0.0001 
Axisymmetric  Modes,  Symmetric,  k  -  0,  h/R  =  0.01 
Axisymmetric  Modes,  Antisymmetric,  k  =  0,  h/R  =  0.01 
n  =  2  Modes,  Symmetric,  k  =  0,  h/R  =  0.01 
n  =  2  Modes,  Antisymmetric,  k  =  0,  h/R  =  0.01 
Axisymmetric  Mode  Shapes,  Symmetric,  k  =  0,  h/R  =  0.01 
Axisymmetric  Mode  Shapes,  Antisymmetric,  k  =  0,  h/R  =  0.01 
Axisymmetric  Modes,  Symmetric,  k  =  0,  h/R  =  0.001 
Axisymmetric  Modes,  Antisymmetric,  k  =  0,  h/R  =  0. 001 
n  =  2  Modes,  Symmetric,  k  =  0,  h/R  =  0.001 
n  =  2  Modes,  Antisymmetric,  k  =  0,  h/R  =  0. 001 
Axisymmetric  Mode  Shapes,  k  =  0,  h/R  =  0.001,  e  =  0. 15 
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NOMENCLATURE 


am’  amplitudes  of  ring  flexural  vibrations 


h 

k 

P 

q 


a 


w 


C,  N,  Q,  S 
Tn,  T, ,  T0  .  .  .  T, 


O’  1’  2 


12 


D 

E 


E  ,  E  ,  E 
a  e  ae 


M  ,  M  .  M 
a’  e  ae 


N  ,  N  .  N 
a  e  ae 


Q  .  Q 


R 


S  ,  S 
a  e 


a 


thickness  of  shell 

(1  -  V2)  K 

pressure 

DQyR2 

(l-e  cos  a)/e 

meridional  displacement  (Figure  1) 
circumferential  displacement  (Figure  1) 
normal  displacement  (Figure  1) 

defined  functions 

Eh3/ 12  (1  -v2) 

Young*  s  modulus 
membrane  strains 

stress  couples 

stress  resultants 

transverse  shear  stress  resultants 

radius  of  the  generating  circle  of  torus  (Figure  1) 
membrane  prestress  forces 

meridional  position  angle  (Figure  1) 
ratio  of  the  two  radii  of  the  torus  (Figure  1) 
viii 


Ka’  *V  Kae 


^a’  ^ ae 


co 

r 

A 


circumferential  position  angle  (Figure  1) 
pR/Eh,  prestress  parameter 
bending  strains 
2  2  2 

( p  R  /E  e  )  to  ,  frequency  parameter 
Poisson' s  ratio 
material  density 

rotations  of  the  normal  to  the  shell 

circular  frequency 
(h/R)2/ 12 

spacing  between  finite  difference  stations 


Matrices 

a,  b,  c 
x 

A,  B,  C,  D,  E,  F,  G,  H,  P 
A,  B,  C 
Z 
Z 


5x5  matrices 
1x5  column  matrix 
4x4  matrices 
3x3  matrices 
1x4  column  matrix 
1x3  column  matrix 


Indices 

i 

I 

m 

n 


station 
last  station 

Fourier  index  of  ring  flexural  vibrations 
Fourier  index  of  shell  vibrations 
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Section  1 
INTRODUCTION 


A  study  of  the  free  vibrations  of  the  pressure  prestressed  toroidal  mem¬ 
brane  [1]  shows  that  the  frequencies  can  be  grouped  into  four  families:  a 
family  of  flexural  vibrations  associated  with  relatively  low  frequencies, 
and  three  families  of  vibrations  with  high  frequencies.  The  flexural  vi¬ 
brations  are  found  to  depend  strongly  on  prestress.  The  vibrations  of  the 
other  three  families  are  practically  insensitive  to  prestress.  Since  both 
the  effects  of  prestress  and  the  effects  of  wall  bending  stiffness  enter  the 
fundamental  equations  in  a  similar  manner  (through  terms  whose  coeffi¬ 
cients  are  small  compared  to  the  coefficients  associated  with  the  membrane 
terms),  it  can  be  expected  that  bending  stiffness  has  a  strong  influence  on 
the  frequencies  of  the  flexural  vibrations.  This  report  presents  the  results 
of  a  numerical  study  of  the  effect  of  wall  bending  stiffness  on  the  frequencies 
and  mode  shapes  of  the  flexural  vibrations. 

Several  iterative  numerical  methods  suitable  for  the  analysis  of  the  free 
vibrations  of  shells  of  revolution  have  appeared  recently  in  the  literature 
[1]  [2]  [3],  Kalnins  [2]  uses  a  multisegment  direct  numerical  integration 
approach  and  Reference  [1]  uses  finite  differences  to  evaluate  a  certain 
determinant  corresponding  to  a  trial  value  of  the  frequency.  Cohen  [3]  uses 
a  method  which  iterates  on  the  mode  shape  instead  of  the  frequency.  Al¬ 
though  the  methods  of  References  [2]  and  [3]  are  applied  to  shells  of  revolu¬ 
tion  without  prestress,  their  extension  to  shells  with  prestress  appears  to 
be  straightforward.  None  of  these  methods,  however,  possesses  a  sub¬ 
stantial  advantage  over  the  others.  This  study  uses  the  numerical  method 
of  Reference  [1], 


Section  2 

FUNDAMENTAL  EQUATIONS 


The  present  analysis  of  the  free  vibrations  of  the  pressure  prestressed 
toroidal  shell  is  based  upon  a  set  of  equations  which  result  from  the 
addition  of  the  bending  terms  of  the  Sanders  linear  shell  theory  [4]  to  the 
prestressed  membrane  equations  presented  in  Reference  [1].  These  are 
(refer  to  Figure  1): 

Equilibrium 
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Constitutive  Relations 
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The  above  constitute  twenty  equations  for  the  twenty  unknowns  N  ,  N  , 

Ot  u 


^of’  R ae ’  ^a’  ^0’  Ka’  Ke’  Kae’  ^ a ’  ^e’  ^a’  ^0’  ^ ae , 


u,  v,  and  w.  The  stress  resultants  S  and  are  known  functions  of 

a  0 


prestress. 


The  state  of  prestress  is  determined  from  a  separate  analysis  of  the 
toroidal  shell  subjected  to  static  internal  pressure.  An  analysis  based 
on  the  linear  membrane  theory  [5]  gives 

1  -  4  e  cos  a 

Sa  =  pR  1  -  e  cos  a  ( 

se  -  |pR 

Analyses  of  the  toroidal  membrane  under  internal  pressure  based  upon 
nonlinear  theories  [6,  7]  show  that  the  linear  meridional  stress  resultant, 

S^,  is  in  error  by  a  negligible  amount,  but  that  the  linear  circumferential 
stress  resultant,  S  ,  can  be  in  error  by  18%.  The  significant  difference, 
however,  between  the  linear  and  nonlinear  stress  distributions  is  confined 
to  a  small  area  of  the  torus  at  a  =  ±  n/2.  Furthermore,  an  analysis  [8] 
of  the  toroidal  shell  under  internal  pressure  based  upon  equations  which 
consider  wall  bending  stiffness  and  nonlinear  behavior  shows  that  the  stress 
resultants  do  not  differ  substantially  from  those  obtained  from  nonlinear 
membrane  theory.  Hence,  it  appears  that  the  state  of  prestress  determined 
according  to  the  linear  membrane  theory  is  adequately  accurate  for  the 
present  analysis  of  vibrations. 
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Section  3 

REDUCTION  TO  SECOND  ORDER  DIFFERENTIAL  EQUATIONS 


The  solution  of  the  fundamental  equations  (1-20)  is  started  by  separating 
the  variables.  Set 


Na'  *V  Mtr  V  V  u 

Ne’  Es’  Me’  VQo'w 


N  ,  E  ,  M  ,  k  ,  <j>  ,  Ru 
an  an  an  an  van  n 

N  ,  E^  ,  ,  k .  ,  ,Rw 

en  en  0n*  0n  R  n 


cos  n  e 
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Ma0’  W  <V  v 


^aen’  ^aan*  ^en’  ^ae n 
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sin  n  e 


Next,  define 

S  =  ( sin  a)  /r 

C  =  ( cos  a )  /  r 

N  =  n/r 

Then,  use  of  equations  (22)  and  (23)  in  equations  (1-14)  yields 
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2  2 
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(31) 
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(32) 
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(33) 
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^ a 

=  -  w'  +  u 

(35) 

*e 

=  Nw  -  Cv 

(36) 

2U 

=  Nu  +  v'  +  Sv 

(37) 

In  equations  (24-37)  and  subsequent  expressions  the  subscript  n  has  been 
dropped.  Prime  indicates  differentiation  with  respect  to  a . 

At  this  point  the  problem  has  been  reduced  to  the  simultaneous  solution  of 
the  ordinary  differential  equations  (24-37)  and  (15-20).  Our  goal  is  to 
derive  four  simultaneous  second  order  differential  equations  for  u,  v,  <j>a> 
and  q  .  First  eliminate  Q  from  equations  (25)  and  (26)  using  equation  (28). 

Oi  0 

Then,  substitute  the  strain  displacement  relations  (29-34),  two  of  the  rota¬ 
tion  expressions  (36,  37),  and  the  constitutive  relations  (15-20)  into  the 
equilibrium  equations  (24-27).  The  resulting  equations  together  with  (35) 
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constitute  five  equations  for  u,  v,  w,  $  ,  and  qQ.  They  may  be  written  in 
matrix  form  as 

ax"  +  bx1  +  c  x  =  0  (38a,  b,  c,  d,  e) 


where 


The  elements  of  the  a,  b,  c  matrices  are  given  in  Appendix  A. 


We  now  manipulate  equations  (38)  so  that  w  can  be  eliminated  using  equation 
(38c)  (equilibrium  of  forces  in  the  normal  direction).  In  order  that  no  deri¬ 
vatives  higher  than  the  second  derivative  appear  in  the  final  equations, 
eliminate  v"  from  equation  (38c)  using  equation  (38b)  (equilibrium  of  forces 
in  the  circumferential  directions)  as  follows: 


v"  =  (b21  u'  +  c2!  u  +  b22  V'  +  C22  v  +  a23  W"  +  b23  w’  +  C23  W 


b0.  <t>  *  +  C0  .  <b  ) 

24  24  ycr 


Then  equation  (38c)  reads: 


b31u’  +  c31u+b32v'+  c32v  +  a33w"+b33w’  +  c33w+b34V 
+  a34  V  r  V  rSqa=  0 


where  the  coefficients  are  given  in  the  appendix.  From  equation  (35) 


w1  =  u  -  6 

w"  =  u1  -  <j>^  (41) 

After  eliminating  the  derivatives  of  w  using  (41),  equation  (40)  can  be 
written  as 

Qw  =  T.u'  +  T„u  +  T0v'  +  T,v  +  Tc<£  '  +  T ,<j)  -  Tq  *  -  r  Sq  (42) 

^  1  2  3  4  5  a  6va  ~ a  '  ' 
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and  the  derivative  of  w  can  be  written  as 


Qw'  =  TjU"  +  (^ '  +  T2)  u'+  T2‘u  +  T^v"  +  (T3‘+  Tj)  v'+  T4'  v 

+  X4"+  (T'  +  T,)  V  +T'<S>  -rq  rSq  r(C-S2)  q  -  Q'w  (43) 
5  a  5  6 a  6  v  ’  ^  v  ’ 

The  T' s  and  Q  are  defined  in  the  appendix. 

Now,  multiply  equations  (38a,  b,  d,  e)  by  Q  and  eliminate  Qw'  from  them 
with  equation  (43).  Then  multiply  the  resulting  equations  by  Q  and  elimi¬ 
nate  Qw  with  equation  (42).  The  resulting  four  equations  for  u,  v,  <j>a>  and 
q^may  be  written  in  matrix  form  as 

AZ"  +  BZ'  +  CZ  =  0  (44) 


where 


fu  1 


v  j 


and  the  elements  of  A,  B,  and  C  are  given  in  Appendix  A. 


The  basic  equations  (44)  simplify  for  the  case  of  axisymmetric  vibrations, 
n  =  0,  v  =  0.  The  equation  of  equilibrium  of  forces  in  the  circumferential 
direction  is  then  satisfied  identically,  and  equations  (44)  become 

AZ”  +  BZ'  +  CZ  =  0  (45) 


where 


Z  =  4 


l:.\ 


The  elements  of  A,  B,  and  C  are  the  remaining  elements  of  A,  B,  and  C  if 
the  second  row  and  column  are  deleted,  and  N  is  set  equal  to  zero. 
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The  fundamental  equations  have  now  been  reduced  to  the  solution  of  two 
sets  of  second  order  differential  equations  as  follows:  1)  four  equations 
(44)  for  the  general  case,  2)  three  equations  (45)  for  the  axisymmetric 
case.  These  differential  equations  will  be  integrated  numerically  by  a 
method  described  in  Reference  [1].  For  completeness,  the  description 
of  the  numerical  method  is  reproduced  in  the  next  section. 

At  this  point  we  remark  that  the  fundamental  equations  (1-20)  could  be 
more  easily  reduced  to  four  second  order  differential  equations  for 
u,  v,w,  and  in  the  manner  of  Budiansky  and  Radkowski  [9].  However, 
application  of  the  numerical  method  described  in  the  next  section  resulted 
in  a  convergence  of  the  finite  differences  which  was  too  slow  for  practical 
computation. 

The  fundamental  equations  can  also  be  reduced  to  four  equations  for  v,w, 

N  ,  and  M  .  First,  derive  a  set  of  five  equations  for  u,  v,w,N  and  M  , 
at  at  M  ’  a  a 

then  solve  the  equation  of  equilibrium  of  forces  in  the  meridional  direc¬ 
tion  for  u,  and  finally  eliminate  u  from  the  remaining  four  equations.  The 
application  of  the  present  numerical  analysis  to  these  equations  resulted 
in  efficient  and  accurate  results  except  in  the  approximate  range 

2  2 

n  <  \  <  n _ 

2  (1  +  v)  (1  +  e)2  2  (1  +  v)  (1  -  e)2 

In  this  range  the  coefficient  of  u  in  the  equation  of  equilibrium  of  forces 
in  the  meridional  direction  has  a  zero. 
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Section  4 

NUMERICAL  ANALYSIS 


Since  the  geometry  and  prestress  are  symmetrical  about  a-  0  and  a  =  ^ r, 
we  need  to  consider  only  one -half  of  the  torus  corresponding  to  the  range 
0  <  a  <  7r .  Let  this  range  be  subdivided  by  I  +  1  equally  spaced  stations. 
Then  the  spacing  between  stations  is 

A  =  7T  / 1 

and  the  position  angle  for  the  i^  station  is 


a.  = 
i 

th 

The  derivatives  of  Z  at  the  i  n 
difference  formulas 


Z.'  = 

1 

1 

2A 

(Z.  .  - 

x  1+  1 

zi-l> 

Z."  = 

1 

1 

A2 

(Z.  ,  - 

1+  1 

2zi +  zi  - 1> 

(46) 

i  =  0,  1,  2, . I 

station  are  approximated  by  the  central 


With  these  formulas  we  obtain  from  equation  (44)  the  set  of  difference 
equations 


where 


D.  Z.  .  +  E.  Z.  +  F.  Z.  =  0 
li+l  ii  li-l 

i  =  0,  1,  2, 


D.  = 

2 

A.  +  B. 

l 

A 

l  l 

E.  = 

. 

4  A.  +  2  AC. 

l 

A  l  l 

F.  = 

2 

A.  -  B. 

l 

A 

l  l 

(47) 


(48) 
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The  difference  equations  (47)  are  augmented  by  conditions  at  i  =  0  and 
i  =  I.  These  end  conditions  are  obtained  from  considerations  of  conti¬ 
nuity  of  the  displacement  functions,  u,  v,  and  w.  Now,  the  solution  of 
the  fundamental  equations  (1-20),  and  the  reduced  equations  (44)  are  such 
that  either  v  and  w  are  even,  u  is  odd  or  v  and  w  are  odd,  and  u  is  an 
even  function  of  a  with  respect  to  a  =  0  and  n  .  These  two  groups  of  solu¬ 
tions  will  be  denoted  as  symmetric  and  antisymmetric  modes  respectively. 
Hence,  the  end  conditions  are 


Z  ^  —  +  GZj, 

ZI+1  -  GZI-1 


(49) 


where  G  is  the  diagonal  matrix 

’-1 

G  -  *-i 

-i 

and  plus  and  minus  signs  refer  to  symmetric  and  antisymmetric  modes 
respectively. 


The  difference  equations  (47)  together  with  the  end  conditions  (49)  make  up 
a  set  of  homogeneous  algebraic  equations.  The  eigenvalues  of  these  equa¬ 
tions  are  obtained  by  trial  and  error,  using  a  special  Gaussian  elimination 
technique  devised  by  Potters  [10]. 


The  equations  for  this  procedure  are  obtained  as  follows.  Let 


Z. 

l 


z.  . 
1+  1 


(50) 


Substitute  this  expression  into  the  difference  equations  (47),  and  by  com¬ 
paring  the  result  with  (50)  obtain  the  recurrence  relation 


P. 

i 


F.P.  .] 
l  l-l 


1  -  1 
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i  =  1,  2,  3, 


(51) 


Now,  write  the  difference  equations  (47)  at  i  =  0  and  eliminate  Z_^  using 
the  first  of  the  end  conditions  (49).  Again,  by  comparing  the  result  with 
equation  (50)  obtain 

P0  ■  V1  [°o  ±  Fo  °1  <52> 


where  plus  and  minus  signs  refer  to  symmetric  and  antisymmetric  modes 
respectively.  Equations  (52)  together  with  the  recurrence  relation  (51) 
provide  all  the  P' s  up  to  Finally,  write  the  difference  equations 

(47)  at  i  =  1.  Eliminate  Z^  +  ^  using  the  second  of  the  end  conditions  (49), 
and  Zj  _  ^  using  equation  (50).  The  result  is 

[Ej  -  (Fj+  DjG)  PI_1]  Zj  =  HZj  =  0  (53) 


Since  Zj  f  0,  we  must  require  that  the  determinant 


V 


Ej  -  (F  +  Dj  G)  P 


I-  1 


(54) 


vanish.  Equation  (54)  is  effectively  a  frequency  equation.  A  value  of  X 
which  gives  v  =  0  is  a  natural  frequency  of  vibration. 


The  mode  shape  corresponding  to  a  natural  frequency  can  be  calculated 
once  a  Xfor  which  V  =  0  is  obtained.  Set  the  amplitude  of  one  of  the  dis¬ 
placements  at  i  =  I  equal  to  unity  and  calculate  the  remaining  unknowns 
at  i  =  I  from  three  of  the  equations  (53).  Thus,  for  the  symmetric  modes 


fol 


For  antisymmetric  modes 


Z 


I 


1 


(55) 


(56) 
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where 


*al 

H33  H34 

-1 

4 

H31 

H.q  Haa 
43  44  _ 

[H4lj 

The  remaining  Z*  s  can  then  be  calculated  from  equations  (50)  in  the 
reverse  order. 


The  w  displacement  is  obtained  from  equation  (42).  In  finite  difference 
form  the  equation  reads 


w. 

i 


rT  1  oT  tTi  •  ■  i-  Zi  •  ,)  +  Tq  .  (Z0  .  ,  -  Z„  .  ,) 

Q.  [  2A  1  l,i  l,i+l  1,1-1  3,i  2,1  +  1  2,i-l 

5.1  3,1  +  1  3,i-l  4,i  +  l  4,i-l  1 

+  Tr,  .  Z.  .  +  T .  .  Z„  .  +  T/  Z„  .  -  r  S.  Z  .  f 

2.1  l,i  4,1  2,i  6  3,1  l  4, l J 


i  =  l,  2,  3, 


1-1 


For  symmetrical  mode  shapes 


1 


W0  Qn  t  A  (T1,0  Zl,l  +  T5,0  Z3,l  “  r  Z4,l)  +  T4,0  Z2,0^ 


<0 


Wl  Qi  Z1,I-1+  T5,I  z3 , 1  - 1  '  rZ4,I-l)  +  T4,0^ 

and  for  anti  symmetrical  mode  shapes 


w0  =  wi  =  0 


The  procedure  may  be  summarized  as  follows: 

1.  Assume  a  value  of  \  ; 

2.  Calculate  the  elements  of  the  A,  B,  C,  D,  E,  F,  and  P 
matrices  at  all  stations  from  equations  (44),  (48),  (51) 
and  (52); 

3.  Calculate  the  determinant  V  from  equation  (54); 


4. 


Repeat  steps  1-3  and  plot  V  versus  x.  From  this  plot 
determine  a  X  for  which  V  =  0.  This  x  is  a  natural 
frequency; 

5.  Calculate  the  mode  shape  corresponding  to  a  natural 
frequency  from  equations  (50)  and  (55-60). 

The  equations  in  this  procedure  were  programmed  for  the  IBM  7094 
computer. 

A  number  of  natural  frequencies  were  calculated  with  1  =  50  and  100  and 
compared.  The  difference  was  found  to  be  less  than  one  percent.  There¬ 
fore,  all  results  were  obtained  with  I  =  50.  For  this  number  of  finite 
difference  spacings  the  IBM  7094  required  approximately  1.5  seconds  to 
evaluate  one  determinant  V  of  the  general  vibrations. 
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Section  5 
RESULTS 


Results  of  the  present  numerical  analysis  are  contained  in  Figures  2-20. 
These  results  show  the  effect  of  bending  stiffness  on  the  vibration  of  shells 
under  1)  high  prestress,  k  =  0.002,  2)  low  prestress,  k  =  0.0001,  and 
3)  no  prestress,  k  =  0.  All  results  are  for  Poisson' s  ratio  v=  0.  3. 


The  modes  of  the  flexural  vibrations  of  the  prestressed  toroidal  shell  may 
be  thought  to  consist  of  two  types  of  vibrations:  1)  the  modes  of  a  torus 
whose  meridional  curve  (cross-section)  is  not  allowed  to  distort  (overall 
ring  vibrations),  2)  the  modes  of  a  pressure  prestressed  circular  ring 
with  radius  R  (cross-sectional  or  ring  flexural  vibrations).  The  first  type 
is  approximated  by  the  vibrations  of  a  thin  ring  of  radius  R/e  without 
prestress  [11]: 


Bending  modes 

symmetric  (in  plane  bending) 


o  2,2  .  J2 

>  _  2  n  (n  -  1) 

A-  “  €  n 

2  (n  +  1) 

antisymmetric  (out  of  plane  bending) 


x  =  e 


o  2/2  .  .2 

2  n  (n  -  1) 

2  (n^  +  1  +  v) 


Extensional  modes  (symmetric) 


X  =  n  +  1 


Torsional  modes  (antisymmetric) 

1  2 

x  l/i  n  . 

2  ^  +  1  +  v 


(61) 


(62) 


(63) 


(64) 
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The  frequencies  of  flexural  vibrations  of  a  pressure  prestressed  ring  of 
radius  R  are  derived  in  Appendix  B: 

X  =  -y  n?2(2m-2"-1-)  U  +  r  (m2-  1)]  (65) 

e  m  +  1 

The  effect  of  bending  stiffness  on  the  frequencies  of  a  shell  with  high  pre¬ 
stress  are  shown  in  Figures  2-7.  These  results  are  for  n  =  0,  1,  and  2, 
k  =  0. 002  and  ^/R  =  0.  01.  It  can  be  seen  that  the  frequencies  of  the  first 
mode  of  each  type  of  vibration  shown  are  practically  unaffected  by  bending. 
However,  the  effect  of  bending  increases  with  the  mode  number  so  that 
the  increase  in  frequency  of  the  fourth  mode  is  of  the  order  of  10%.  This 

trend  is  forecast  by  the  frequencies  of  the  uncoupled  ring  flexural  vibra- 

2 

tions  given  by  equation  (65)  in  which  the  r  part  increases  with  (m  -  1). 

The  effect  of  bending  is  larger  on  shells  with  high  e  (fat  toroids).  This  is 
also  to  be  expected  because  for  high  e  the  mode  shapes  consist  mostly  of 
local  deformation  of  the  meridional  curve  [1].  For  e  =  0.  75,  n  =  0  sym¬ 
metric  vibrations,  the  mode  shapes  of  the  first  four  modes  of  a  shell  with 
and  without  bending  are  compared  in  Figure  8.  This  figure  shows  that  the 
effect  of  bending  on  these  mode  shapes  is  small.  In  summary,  we  can  say 
that  for  shells  with  k  =  0.  002  and  =  0.  01  the  effect  of  bending  on  the 
first  four  modes  of  the  n  =  0,  1,  2  vibrations  is  only  moderate. 

The  effect  of  bending  on  the  axisymmetric,  symmetric,  frequencies  of  a 
shell  with  k  =  0.0001  and  =  0.01  is  shown  in  Figure  9.  In  this  case, 
bending  has  increased  the  frequencies  by  approximately  a  factor  of  two 
when  6=0.  75.  For  smaller  e  the  increase  is  smaller.  For  a  shell  with 
k  =  0.  0001  and  =  0.  001  the  frequencies  increase  by  less  than  5%. 

Frequency  curves  for  a  shell  without  prestress,  k  -  0,  n  =  0  and  2  and 
h/R  =  0.01  are  shown  in  Figures  10-13.  The  dashed  curves  in  these  fig¬ 
ures  represent  the  frequencies  of  the  uncoupled  modes.  The  shell  frequen¬ 
cy  curves  are  the  solid  curves.  These  curves  exhibit  transition  regions 
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which  occur  as  sharp  jumps  for  the  axisymmetric,  symmetric,  frequency 
curve,  and  as  more  gradual  changes  for  the  axisymmetric,  antisymmetric, 
and  n  =  2  frequency  curves.  In  a  qualitative  sense  these  frequency  curves 
exhibit  the  same  features  as  those  of  the  prestressed  membrane  when 
k  =  0.002  [1]  (or  refer  to  the  dashed  curves  of  Figures  2-7  of  this  report). 

The  axisymmetric  mode  shapes  for  a  shell  without  prestress,  k  =  0,  and 
h/R  =  0.  01  are  shown  in  Figures  14  and  15.  The  first  column  of  mode 
shapes  in  these  figures  shows  essentially  the  uncoupled  mode  shapes  for  e 
below  the  transition  region.  In  this  region  the  overall  ring  mode  shape 
prevails  for  the  symmetric  vibrations,  but  for  the  antisymmetric  vibra¬ 
tions  an  irregular  shape  appears.  The  third  column  again  shows  essen¬ 
tially  the  uncoupled  mode  shapes.  These  are  for  e  above  the  transition 
region  and  are  one  order  of  complexity  higher  than  those  below  the  transi¬ 
tion  region.  The  fourth  column  shows  mode  shapes  for  e  =  0.  75  which  are 
predominantly  local  oscillations  of  the  meridional  curve.  In  a  qualitative 
sense,  the  mode  shapes  of  the  shell  without  prestress  and  ^/r  =  0.  01 
exhibit  the  same  features  as  those  of  the  prestressed  membrane  when 
k  =  0.002  [1]. 

Frequency  curves  for  a  shell  without  prestress,  k  =  0,  n  =  0  and  2  and 
h/R  =  0.001  are  shown  in  Figures  16-19.  The  significant  feature  of  these 
curves  is  their  tendency  to  pair  up  with  nearly  the  same  frequencies. 
Examples  of  the  mode  shapes,  shown  in  Figure  20,  associated  with  fre¬ 
quencies  that  pair  up  appear  to  be  mirror  images  of  each  other.  That  is, 
if  the  mode  shape  of  the  second  symmetric  mode  is  folded  on  top  of  the 
mode  shape  of  the  first  symmetric  mode,  then  the  two  mode  shapes  are 
nearly  the  same.  The  third  and  fourth  symmetric,  second  and  third  anti¬ 
symmetric  and  fourth  and  fifth  antisymmetric  mode  shapes  shown  in  Fig¬ 
ure  20  are  mirror  images  also.  Another  significant  feature  of  the  mode 
shapes  presented  in  Figure  20,  is  the  fact  that  the  motion  takes  place 
mainly  near  the  crowns. 
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The  paired  up  frequencies  with  nearly  the  same  mode  shapes  can  be  ex¬ 
plained  by  the  following.  Let  the  radius  of  the  torus  go  to  infinity  (e -►()). 
Then  the  torus  approaches  a  cylindrical  shell.  At  every  natural  frequency 
of  the  cylindrical  shell  there  are  actually  two  identical  frequencies  with 
identical  mode  shapes.  These  frequencies  separate  slightly  and  the  mode 
shapes  become  slightly  different  when  imperfections  (for  example,  when 
the  cylindrical  shell  is  not  exactly  circular)  are  present.  Thus  the  pairing 
up  of  the  frequencies  of  the  toroidal  shell  are  due  to  the  weak  influence  of 
the  circumferential  curvature  which  has  the  equivalent  effect  of  imperfections 
on  the  frequencies  of  the  cylindrical  shell. 
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Section  6 
CONCLUSIONS 


The  following  conclusions  are  drawn  from  this  study: 


1)  Bending  stiffness  has  only  a  moderate  effect  on  the  natural 
frequencies  and  mode  shapes  of  relatively  thick  shells, 
(%  =  0.01)  under  high  prestress  (k  -  0.002). 

2)  When  prestress  is  small  (  «■  -  0.0001)  bending  stiffness 
increases  the  natural  frequencies  of  relatively  thick  and 
fat  toroidal  shells  (e  large,  =  0.01)  by  as  much  as  a 
factor  of  two. 


3)  In  the  absence  of  prestress,  when  membrane  theory  pre¬ 
dicts  a  continuous  frequency  spectrum  with  discontinuous 
mode  shapes,  the  consideration  of  bending  stiffness  leads 
to  discrete  frequencies  with  continuous  mode  shapes. 


4)  Bending  stiffness  should  be  considered  in  the  calculation 
of  the  natural  frequencies  and  mode  shapes  associated 
with  the  flexural  vibrations  of  toroidal  shells.  ^ 


{A 


v  - 
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APPENDIX  A 


Elements  of  Matrices 

The  non  zero  elements  of  the  a,  b,  and  c  matrices  are: 

an  =  l  +  k(l+|c) 

a22  =  |(l-v)  +  k(l+ic)  +  i  (1-v)  r  (1  +  3C)2 

a23  =  'I(1'v)  r<1  +  3C>N  =  -  a32 

a33  "  2  ^  rN 

a54  =  1 

bu  =  (l+jk)S 

b12  =  N[|(l-v) -|(l-v)  r(i+C)(l  +  3C)  ]  =  -b21 

b13  =  1  -  vC  +  k(l+|c)  +  ia-v)  r(l  +  C)N2 

b22  =  |S[  l-v  +  k-i(l-v)  r  (1  +  3C)  (5  +  3C)  ] 

b23  =  |(1'v)  r  N  S  (2  +  3  C) 

b24  =  r  N  [  vC+|(l-v)  (1+3C)  ]  =-  rb52 

b31  =  1  -  vC  +  ±kC  +  i(l  -  v)  r(l+C)  N2 

b32  =  ’  |(1  -  v)  rNS 

b33  =  j  (l  -  v)  r  n2S 

b34  =  k  (1+ic) +  i(l  +  v)  rN2 
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b35 

— 

-  r 

b43 

= 

1 

b53 

= 

|(l  +  v)  N2 

b54 

= 

S 

C11 

= 

(l-v2)e2\+  vC-S2-i  (l-v)N2-ik(N2  +  S2)  -  i  (1  -  v)r(l  +  C)2N2 

C12 

= 

NS  [-|(3-  v)  -  k  +  |(l-v)r(l+C)  (1  +  3C)] 

c13 

= 

S  (l  +  C)  [l+|k-|(l-v)rN2] 

C14 

= 

-|[kC  +  ±(l-v)r(l+C)  N2] 

C15 

= 

r 

C21 

= 

-NS[|(3-v)  +  k  +  |(l-v)r(i+C)2] 

C22 

= 

(1  -  v2)  e2x  -  |(1  -  v)  (C  +  S2)  -  N2  -  ^  k  (N2+  S2+  C2)  -  kC-  rN2C2 

+  ~  (1  -  v)  r  (1  +  3  C)  [S2  (5  +  3  C)  -  C  (1  +  3  C)] 

C23 

= 

N  |  -  v  +  C  +  k(l+C)  +  rN2C  +  ~(l-v)r[C(l  +  3C)  -  3S2(l+C)]) 

C24 

= 

r  N  S  [C  - 1  (1  -  v)] 

C31 

= 

S  [v-C-  k  (1+jC)] 

C32 

= 

-  N  |  -  v  +  C+  k  (l  +  C)  +  TN2C  -  i(l  -  v)  r  [S2-  |c  (1  +  3C)]J 

C33 

= 

-  (1  -  V2)  e2x  +  1  -  2vC  +  C2+  jkC  (3  +  C)  +  jkN2+  rN4+  (1  -  v)TN2C 

C34 

= 

|s  [k  +  (3-  v)  rN2] 

C35 

= 

-  r  S 

22 


C51  =  -^(l-v)(l+C)N2 

c52  =  NS  [C  +  v  (1+C)  +^(1  -  v)  (1  +  3 C)] 

c53  =  '2n2s 

c54  =  V  C  -  S2  -  \  (1  -  v)  N2 


The  coefficients  in  equation  (40)  are: 


a33  a33  "  T0  a23 


b3i  b3i  '  T0  b2i 


C3i  C3i  "  T0  C2i 


i  =  1,  2,  3,  4 
i  =  1,  2,  3,  4 


where 


T  = 

r\ 


i(l-v)r(l  +  3C)N 


0  |  (l-v)  +  k(l+ic)  +  i(l-v)r(l  +  3C)i 


The  elements  of  the  A,  B,  and  C  matrices  are: 
An  =  Q2  [1  +  k  (1+^C)]  +  QTL  T? 

A12  =  ^  T3  T7 

A13  =  Q  T5  T7 

A14  =  'rQT7 
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Bn  =  Q2S(l+ik)  +  QT7(T1'  +  T2)  +  T1  Tg 
.  B12  =  |Q2N  [(1  +  v)  -  \  (1  -  v)  r  (1+C)  (1  +3C)]  +  QT?  (Tg?  +T4)  +  Tg  Tg 

B13  =  QT7(T5,  +  V  +  T5T8 

B14  =  -  r  (Q  s  T?  +  Tg) 

B2i  =  -|Q2N[l  +  v+i(l-v)r(l  +  3C)(l-C)]  +  QT9(T1'  +  T2)  +  T1T10 
B22  =  ^Q2s  [1  -  V  -  \  (1  -  v)  r  (1  +  3C)  (5 +  3C)]  +  Q  T9  (Tg-  +  T4)  +  Tg  TJ0 

B23  =  rQ2N[vC+i(l-v)(l+3C)]  +  QT9(T5'  +  T6)  +  T5T10 


24- 


^24  _  "  r  (Q  S  +  T10) 

B31  =  Q  <Ti'  +  T2>  -  Q'  T1 

B32  =  Q(T3'  +  T4)-Q<T3 

B33  =  Q  (V  +  T6>  '  Q*  T5 

B34  =  -  r(Qs  '  Q') 

B41  =  Q  T11  +  T2)  +  T1  T12 

B42  =  -  Q2n  [vC  +{(1  -  v)  (1+3C)]  +  Q  Tn  (T3-  +  T4)  +  T3  Ti2 

b43  -  Q2s  +  QTU<TS'*T6)  +  T5T12 

B44  =  -r(QST^+T^2) 

Cu  -  Q2[(l-v2)e2x  +  vC-S2-i(l-v)N2-ik(N2+S2)-i(l-v)r(l+C)2N2] 

+  Q  t7  t2'  +  t2  t8 

C12  =  Q2NS  [  -  j(3  -  v)  -  k  +  g  (1  -  v)  r  (1+C)  (1+3C)  +  QT?  T4'  +  T4  T8 

C13  =  -jQ2[kC+i(l-v)r(l+C)N2]  +  QT7T6-  +  T6T8 

C14  =  r[Q2-QT7(C-S2)-ST8] 

C21  =  -Q2NS[i(3-v)  +  k+|(l-v)r(l+C)2]  +  QT9T2'  +T2T10 

C22  -  Q2{u-v2>t2x-i(l-v)(C+S2)-N2-|k(N2+S2+C2)-kC-rN2C2 

+  i(l-v)r  (1  +  3C)  [S2  (5  +  3C)  -C(1+3C)]}  +  QT9T4’  +  T4T10 

C23  =  rQ2NS[C-i(l-v)]  +  QT9T6'  +T6T10 
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where 


C24  =  -r  [QT9(C-  S2)  +  ST10] 

Si  =  'Q2  +  Qt2'  -Q't2 

C32  =  QT4'-Q,T4 

C33  =  Q2+QT6'-Q’T6 

c34  =  -  r  [Q  (C  -  S2)  -  Q'  S] 

C41  =  -  jU-v)  (l+C)N2Q2+QTn  T2,+  T2T12 

C42  =  Q2NS[C+  v(l+C)+i(l-v)(l+3C)]  +  QTuT4'+T4T12 

C43  =  Q2[vC-S2-I(l-v)N2]  +  QTnT6«+T6T12 

C44  =  -Q2-r[QTn  (C-S2)  +  ST12] 

Tj  =  1-vC+ikC-  ~  (l-v)rN2(l-C)  +|NT0[l+v+i(l-v)r(i+3  C)(l -C)] 
T2  =  S  [v-C-k(l+|c)  +|(l-v)rN2]  +  NST0[|(3-v)  +  k 
+  (l-v)r  (l  +  C)2  -  i  (l-v)r  (2  +  3C)] 

T3  =  "  rNS  -  |STQ  [1-v  +  k  +  |(l-v)r(l+3C)(5+  3C)] 

T4  =  n{  v-C-k(l+C)  +  r[-N2C+i(l-v)S2-  i(l-v)(l+3C)  C] } 

-  Tq£  (l-v2)e2X'2-(l-v)(C  +  S2)-N2-|k  (N2+S2+C2)  -  kC 
-  r  N2C2+i  (l-v)r  (1  +  3C)  [S2  (5  +  3C)  -C(l  +  3C)]j 
T5  =  k  (1  +  jC)  +  rN2-rNTQ  [vC+|(l-v)  (1+3C)] 

T6  =  s{  |k  +rN2+  ^TNTq  [1-v  +  (1-3  v)  C]J 

Tj  =  1  -  vC  +  k  (1  +  ^C)  +  i  (1-v)  r(i+C)  N2 
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Tg  =  QS(l+C)[l  +  ik-i(l-v)rN2]-Q'T7 
T9  =  i(l-v)TNS  (2  +  3 C) 

T10  =  QN  [-v  +  C  +  k  (1+C) +  TN2C +^(l-v)r(l  +  3C)(C-2S2) 

-i(l-v)r(l-3C)  S2]  -  Q'T9 
Tn  -  i(l  +  v)N2 
T12  =  -  (2N2SQ  + Q'Tn) 

and  prime  indicates  differentiation  with  respect  to  a.  In  differentiating 
N,  S,  and  C  the  following  formulas  are  useful: 

N'  =  -  NS 
S'  =  C  -  s2 
C'  =  -  S  (1+  C) 
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APPENDIX  B 


Flexural  Vibration  of  the  Prestressed  Circular  Ring 


The  vibrations  of  a  circular  ring  of  radius  R,  prestressed  by  pressure  p, 
may  be  determined  from  equations  (15-21)  and  (24-37)  by  setting  e  =  v  =  v 
=  =  0.  This  results  in  the  following  equations: 

Equilibrium 


N'+4q  +  S  (E  •  -  *  )  +  p  R  4  +  p  h  R2  a:2  u  =  0 
a  R  a'  a  ycr  F  K 

N  -4Q'+S  U  •  +  E  )  -  pRE  -  phR2co2w=0 
a  R^a  a'va  a  v  a 

M  '  -  4  Q  =  0 
a  R  xa 

Strain  Displacement 

E  =  u'  +  w 
a 

Rk  —  6  1 
a  va 

6  =  -  w'  +  u 

a 

Constitutive  Relations 


(Bl) 

(B2) 

(B3) 


(B4) 

(B5) 

(B6) 


EhE 

a 

Eh3 
12  Ka 


(B7) 

(B8) 


Prestress 

S  =  pR 
a  v 

Substitution  of  (B3-B9)  into  Bl  and  B2  yields 

2 

(1  +  k+  r)  u"  +  e  x  u  -  rw' 11  +  (1  +  k)  w'  =0 

2 

-  Fu"  '  +  (1  +  k)  u'  +  rw'  "  '  -  *wM  +  (1  -  e  x)w  =  0 


(B9) 
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(BIO) 

(Bll) 


which  with 


u  =  a  sin  m  a 
m 

w  =  b  cos  m  a 
m 

yields  the  following  frequency  equation 

( e2X)2  -  [1  +  m2  +2  k  m2  +  r  m2  (1  +  m2)]  ( e^X) 

+  m2  (m2- 1)  [  k  (1  +  k)  +  Krm2+  r(m2-  1)]  =  0 


For  k  <<1  and  m  not  large  the  lowest  root  of  this  frequency  equation 
is  approximately 

1  m2(m2-l)  r  ^,2 

x  «  ~2  - ^2 - -  U+r(m  "I)] 

t  m  +  1 


(B12) 

(B13) 


(B14) 


(B15) 
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AXIS  OF 
ROTATION 


FIGURE  |:  GEOMETRY  AND  NOTATION 


FIG.  2s  AXISYMMETRIC  MODES ,  SYM METRIC,  K  =  0.002 


FIG.  3:  AXISYMMETRIC  MODES,  ANTISYMMETRIC, /<  =  0.002 


FIG. 4:  n^l  MODES,  SYMMETRIC,  /C  =  0.002 


FIG. 6:  n-2  MODES,  SYMMETRIC,  K=  0.002 


cr  u. 


FIG. 9:  AXISYMMETRlC  MODES,  SYMMETRIC, /C  =  O.OOOI 


FIG.I4:  AXISYMMETRIC  MODE  SHAPES,  SYMMETRIC,  K=0  ,  h/R  =0.01 
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AXIS  OF  ROTATION 


€=O.OOI  €=  0.005  €=0.1  €=0.75 

X=  0.497  X=  0.417  X=  0.0145  X  =  0.00239 


€=0.055  €  =  0.075  €  =  0.15  €  =  0.75 

X=  1.183  X=  1.037  X=  0.753  X=  0.239 


FIG.  15;  AXISYMMETRIC  MODE  SHAPES,  ANTISYMMETRIC,  K-0,  h/R  =  O.OI 
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AXIS  OF  ROTATION 


1-2 


NASA-Langley,  1965  CR-296 


FIG.  20:  AXISYMMETRIC  MODE  SHAPES,  K  =  0  ,  h / R  =0.001,  €  =  0.15 


